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ABSTRACT 

We present measurements of the star formation efficiency (SFE) in 3D numerical simulations of driven 
turbulence in supercritical, ideal-MHD, and non-magnetic regimes, characterized by their mean normal- 
ized mass-to- flux ratio /x, all with 64 Jeans masses and similar rms Mach numbers (~ 10). In most cases, 
the moderately supercritical runs with fi = 2.8 have significantly lower SFEs than the non-magnetic 
cases, being comparable to observational estimates for whole molecular clouds (< 5% over 4 Myr). Also, 
as the mean field is increased, the number of collapsed objects decreases, and the median mass of the 
collapsed objects increases. However, the largest collapsed-object masses systematically occur in the 
weak-field case n = 8.8. The high-density tails of the density histograms in the simulations are de- 
pressed as the mean magnetic field strength is increased. This suggests that the smaller numbers and 
larger masses of the collapsed objects in the magnetic cases may be due to a greater scarcity and lower 
mean densities (implying larger Jeans masses) of the collapse candidates. In this scenario, the effect of a 
weak field is to reduce the probability of a core reaching its thermal Jeans mass, even if it is supercritical. 
We thus suggest that the SFE may be monotonically reduced as the field strength increases from zero to 
subcritical values, rather than there being a discontinuous transition between the sub- and supercritical 
regimes, and that a crucial question to address is whether the turbulence in molecular clouds is driven or 
decaying, with current observational and theoretical evidence favoring (albeit inconclusively) the driven 
regime. 

Subject headings: ISM: clouds — MHD — Stars: formation — Turbulence 



1. INTRODUCTION 

The fraction of a molecular gas mass that is converted 
into stars, known as the star formation efficiency (SFE), 
is known to be small, ranging from a few percent for en- 
tire molecular cloud complexes (e.g., Myers et al. 1986) to 
10-30% for cluster-forming cores (e.g., Lada & Lada 2003), 
even though molecular clouds in general have masses much 
larger than their thermal Jeans masses, and should there- 
fore be undergoing generalized gravitational collapse if no 
other processes prevented it (Zuckerman & Palmer 1974). 
Thus, this reduction of the mass that is deposited in col- 
lapsed objects needs to be accounted for by models of the 
star formation process. 

In the so-called "turbulent" model of star formation 
(see, e.g., the reviews by Vazquez-Semadeni et al. 2000; 
Mac Low & Klessen 2004; Vazquez-Semadeni 2004), the 
low efficiency arises because the supersonic turbulence 
within the clouds, while contributing to global support, 
generates large-amplitude density fluctuations (clumps 
and cores) , some of which may themselves become locally 
gravitationally unstable and collapse in times much shorter 
than the cloud's global free-fall time. Thus, collapse occurs 
locally rather than globally, and involves only a fraction of 
the cloud's total mass. This fraction constitutes the SFE, 
and depends on the global properties of the cloud, such 
as its rms Mach number, the number of Jeans masses it 
contains, and the turbulence driving scale (Leorat, Passot 
& Pouquet 1990; Klessen, Heitsch & Mac Low 2000), or, 
perhaps more physically, on the scale at which the turbu- 
lent velocity dispersion becomes subsonic (Padoan 1995; 



Vazquez-Semadeni, Ballesteros-Paredes & Klessen 2003). 
Nevertheless, the SFE appears to still be too large in non- 
magnetic configurations, being ~ 30% in simulations with 
realistic parameters (e.g., Klessen, Heitsch & Mac Low 
2000; Vazquez-Semadeni, Ballesteros-Paredes & Klessen 
2003), and it is important to investigate the contribution 
of the magnetic field in further reducing the SFE. This 
remains an open issue. In three-dimensional (3D) simu- 
lations of driven, self-gravitating, ideal MHD turbulence, 
Heitsch, Mac Low & Klessen (2001) studied the evolution 
of the mass fraction in collapsed objects (a measure of the 
SFE) in supercritical and non-magnetic cases. They re- 
ported, however, that any systematic trends with the field 
strength that might have been present in their simulations 
were blurred by statistical fluctuations from one realiza- 
tion to another. 

More recently, Li & Nakamura (2004, hereafter LN04) 
and Nakamura & Li (2005, hereafter NL05), have mea- 
sured the SFE in two-dimensional (2D) simulations of de- 
caying turbulence including ambipolar diffusion (AD) and 
(in NL05) a model prescription for outflows. LN04 found 
that the initial turbulence can accelerate the formation 
and collapse of cores within the clouds. NL05 concluded 
that SFEs comparable to those of whole molecular clouds 
(a few percent) required moderately subcritical conditions, 
while moderately supercritical cases gave efficiencies com- 
parable to cluster- forming cores (~ 20%). However, since 
their simulations were done in a decaying regime and in a 
closed numerical box, the SFEs they measured are prob- 
ably upper limits. Real clouds may not be in a decaying 
regime (Vazquez-Semadeni et al. 2005, hereafter Paper I; 
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see also the discussion in §4) and moreover probably un- 
dergo partial dispersion in response to their turbulent en- 
ergy contents, as in the simulations of Clark & Bonnell 
(2004) and Clark et al. (2005), thus reducing the amount 
of mass available for collapse. 

Paper I studied the formation, evolution and collapse of 
the cores formed in 3D simulations of driven MHD tur- 
bulence, albeit neglecting AD. Due to this setup, no col- 
lapse could occur in subcritical cases. Although no quan- 
titative measurements of the SFE were reported there, a 
trend towards decreasing collapse rates with increasing 
field strengths was clearly observed in supercritical and 
non-magnetic cases. This was not due to longer individ- 
ual core collapse times in the magnetic cases with respect 
to the non-magnetic one, since the timescales for formation 
and collapse of the cores were similar in both cases. In- 
stead, the reduction of the SFE was apparently due to a re- 
duced formation rate of collapsing objects in the magnetic 
simulations in comparison with the non-magnetic case. 

The goal of the present paper is to report quantitative 
measurements of the SFE and of the numbers and masses 
of collapsed objects forming in the simulations of Paper 
I and three new sets of similar ones, as a function of the 
mean field strength. In order to overcome the difficulties 
encountered by Heitsch, Mac Low & Klessen (2001), we 
study these variables at fixed settings of the random tur- 
bulence driver. 

2. NUMERICAL METHOD, SIMULATIONS AND PROCEDURE 

We refer the reader to Paper I for details on the simula- 
tions and resolution considerations. We consider four sets 
of simulations at a resolution of 256 3 keeping all physical 
parameters constant, except for the mass-to-magnetic flux 
ratio if, which we vary to investigate the effects of the mag- 
netic field. One set consists of the four simulations pre- 
sented in Paper I with rms Mach number M « 10, Jeans 
number J = L/Lj = 4 (where L is the numerical box size 
and Lj is the Jeans length), and /i = 0.9, 2.8, 8.8 and 
oo, corresponding to subcritical, moderately supercritical, 
strongly supercritical and non-magnetic cases. We refer 
to this set by the label "Paper I" . We also consider three 
more sets of three similar simulations each (with fi = 2.8, 
8.8 and oo), but varying the seed of the random turbulence 
driver. We label the sets by their seed numbers, as "Seed = 
0.1", "Seed = 0.2", and "Seed = 0.3". The driving is com- 
puted in Fourier space, and applied at the largest scales in 
the simulation (<~ 1/2 of the box length), so that it is not 
expected to be the main driver of the local evolution of 
the clumps and cores, because, on the scales of the cores, 
the applied force is nearly uniform, and its main effect 
should just be to push the cores around without severely 
distorting them. 

Although the simulations are scale-free, for reference, a 
convenient set of physical units is n = 500 cm~ 3 , u a = 
c s = 0.2 km s _1 , Lq = L = 4 pc, and to = Lq/uq = 20 
Myr. The latter is the sound crossing time across the box. 
Taking the mean molecular mass as m — 2.4toh, the nu- 
merical box then contains 1.86 x 10 3 M Q . The mean field 
strengths for the /j, = 0.9, 2.8, 8.8 and oo cases are respec- 
tively B = 45.8, 14.5, 4.58, and /xG, corresponding to 
values of [3, the ratio of thermal to magnetic pressure, of 
P = 0.01, 0.1, 1 and oo. The simulations are run for 0.5 



code time units (10 Myr) before turning on the self-gravity. 

As a measure of the SFE, we consider the evolution of 
the collapsed mass fraction M co \ of the simulations. By 
this we mean gas that is at densities n > 500n , since in 
Paper I we noticed that once an object reached densities 
~ 300no it was always already on its way to collapse, a fact 
which is confirmed by the fact that these objects never dis- 
perse during the subsequent evolution of the simulations. 
Throughout the paper, we refer exclusively to collapsed 
objects rather than stars, because our spatial resolution 
is clearly insufficient to determine whether a collapsing 
object eventually breaks up into more fragments to form 
several stars. This is likely to be the case for the most 
massive collapsed objects, with masses up to ~ 100 M . 
Thus, the masses reported in fig. 2 should not be neces- 
sarily interpreted as individual stellar masses, and may 
well be cluster masses. Also note that in runs that form 
a single collapsed object, the plots of M co \ vs. time really 
represent the accretion history onto that object, rather 
than the continuous formation of new collapsed objects. 

As a representative cloud lifetime we take t c \ = 4 
Myr, a time scale that agrees with the estimate of 
Bcrgin ct al. (2004) of 3-5 Myr and also with the 
two-turbulcnt-crossing-times criterion used by Vazquez- 
Semadeni, Ballesteros-Paredes & Klessen (2003). Quoted 
values of the SFE refer to the collapsed mass fraction at 
t = t c i. Note, however, that, because of numerical prob- 
lems when the density contrast becomes too large, not all 
simulations reach this time, although this will not be a 
limitation for the conclusions we will draw. 

3. RESULTS 

Figure 1 shows the evolution of the accreted mass frac- 
tion for the four sets of runs we consider. The moder- 
ately supercritical runs, with /z = 2.8, have generally lower 
SFEs than both the strongly supercritical (/z = 8.8) and 
non-magnetic (p, = oo) runs. In turn, the /i = 8.8 cases 
have SFEs that are generally very similar to those of the 
non-magnetic runs, except in the runs from Paper I. It is 
also noteworthy that the n = 2.8 cases have SFEs « 0.04, 
0.12, 0.025 and 0.05, respectively for each set. Thus, in 
three out of the four statistical realizations, M co \ < 5% at 
t = r c i, in reasonable agreement with the observed SFEs 
at the level of global molecular clouds. 

What is even more interesting is the different way in 
which the magnetic and the non-magnetic runs reach their 
respective collapsed fractions. Figures 5 and 9 in Paper I, 
and their corresponding animations in the electronic ver- 
sion, show that, while the magnetic runs do so with one or 
two relatively massive collapsed objects, the non-magnetic 
run does so with several objects, many of them with low 
masses. The same trend is observed in the three additional 
sets of statistical realizations considered in the present pa- 
per. In fig. 2 we show the masses of the individual col- 
lapsed objects (n > 500no) for all the runs we consider 
at t = 4 Myr, except for those cases in which the simu- 
lation terminated prematurely, in which case we plot the 
collapsed object masses at the last timestep of the simula- 
tion. We see that the non-magnetic runs typically produce 
many more collapsed objects than the magnetic runs, and 
with mass distributions that extend to significantly lower 
values. On half the cases, the minimum masses of the col- 
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lapsed objects increase monotonically with increasing field 
strength. 

One of the main results of Paper 1 was that the core for- 
mation+collapse time scale was not significantly different 
between the magnetic and non-magnetic cases. Thus, the 
above result on the masses and numbers of the collapsed 
objects suggest that the decreased efficiency of the mag- 
netic cases in comparison to the non-magnetic ones arises 
from a decreased probability of collapse events as the field 
strength increases, rather than from an increase of the col- 
lapsing object lifetimes. This suggestion is supported by 
the probability distribution of the density fluctuations for 
the various runs. Figure 3 shows the histograms of the 
density values in each of the simulations, averaged over 
the last five snapshots before turning on gravity. In this 
figure, the panel for the runs from Paper I includes the 
subcritical run with fi = 0.9, in order to see the effect of 
the magnetic field in this case as well, even though this 
run did not undergo collapse. 

Figure 3 shows a clear trend towards decreasing width 
and lower high-density tails with increasing mean field 
strength (increasing fi) , at least over the range of magnetic 
field strengths we have considered here. That is, the prob- 
ability of producing large density enhancements decreases 
with increasing magnetic field strength. This is in agree- 
ment with previous results by Passot, Vazquez-Semadeni 
& Pouquet (1995), Ostriker, Gammie & Stone (1999), 
Hcitsch, Mac Low & Klessen (2001) and Ballcsteros- 
Paredes & Mac Low (2002). Note, however, that in those 
works a reverse trend towards increasing fluctuation ampli- 
tude was observed at stronger values of the field, which we 
do not observe here, again at least over the field strength 
range spanned by our simulations. 

The trend towards lower and less extended high-density 
tails in the histograms at larger field strengths is consistent 
with the trend of the SFE to decrease and of the minimum 
collapsed masses to increase in the same limit, since fewer 
density fluctuations can reach the threshold for collapse, 
and simultaneously, the local values of the Jeans mass are 
larger, so that, in order to collapse, an object needs to 
acquire more mass. In this picture, the field's effect on 
the SFE is only indirect, through its modification of the 
density histogram, rather than by directly increasing the 
minimum mass for collapse (i.e., by causing the magneti- 
cally critical mass of the cores to be larger than their Jeans 
masses). Indeed, in the cases studied in Paper I, examples 
of both collapsing and non-collapsing cores were supercrit- 
ical, and the occurrence of collapse depended on whether 
they acquired the Jeans mass. Similarly, Li ct al. (2004) 
found that all the cores in their supercritical simulations 
were supercritical. 

Finally, it is worth noting that the most massive ob- 
jects systematically arise in the strongly supercritical cases 
(/i = 8.8) in all four simulation sets. The origin of this ef- 
fect, as well as a test of the mechanism suggested above for 
reducing the SFE and increasing the minimum masses will 
require detailed measurements of the field morphology and 
the evolution of the energy balance in the cores prior and 
during the onset of collapse, to be presented elsewhere. 



4. DISCUSSION 



The above results can be placed in the context of previ- 
ous studies. We have found that the presence of a magnetic 
field can further reduce the SFE with respect to the non- 
magnetic case, even in supercritical configurations. Pre- 
vious successful determinations of the effect of the mag- 
netic field have been restricted to decaying, 2D simula- 
tions (LN04; NL05). These authors found that values of 
the SFE comparable to those observed in entire molec- 
ular clouds (a few percent; Myers et al. 1986) required 
subcritical environments and AD-mediated collapse, while 
supercritical environments gave values of the SFE closer 
to those of cluster-forming cores (~ 15% after 1 global 
free-fall time). 

Instead, in our simulations we have found that a mod- 
erately supercritical environment and reasonably realistic 
values of the Mach number (10) and of the Jeans num- 
ber (64 Jeans masses in a 4-pc cube) already give SFEs 
<5% in 3/4 of the cases studied after 0.8 global free-fall 
times (4 Myr), in spite of having a mass-to- flux ratio more 
than twice larger than theirs. The difference is probably 
due mainly to the choice of global setup (2D, decaying, 
versus 3D, driven), since a continuously-driven simulation 
maintains the turbulent support throughout its evolution, 
while a decaying one loses it over time. 

In the simulations of LN04 and NL05, the main role of 
turbulence is to accelerate the initial formation of the cores 
and other density structures, which occurs on the turbu- 
lent crossing time rather than in the long AD timescalc. 
However, as the turbulence decays, its role in providing 
support and producing further fragmentation (which in- 
duces local collapse involving small fractions of the total 
mass) progressively decreases. Indeed, the turbulent Mach 
numbers in the simulations of NL05 had already decayed 
to values ~ 2-3 by the times the collapsed objects were 
forming. Thus, in their supercritical cases, the entire bulk 
mass of the simulation is in principle available for collapse 
at long enough times, although the residual turbulent frag- 
mentation still drives local collapse events first. 

This suggests that a fundamental question in under- 
standing the SFE in molecular clouds is whether real 
molecular clouds are driven or decaying, and sub- or super- 
critical. Concerning the former dichotomy, observational 
evidence tends to suggest that molecular clouds are driven, 
as we discussed at length in Paper I; a few additional con- 
siderations are as follows. If the turbulence is generated 
by instabilities in the compressed layers as the clouds are 
forming (e.g., Hunter et al. 1986; Vishniac 1994; Walder 
& Folini 2000), the injection of turbulent energy is likely 
to last for as long as the accumulation process lasts. Af- 
terwards, the cloud is likely to disperse, as indicated by 
the facts that the gas has disappeared from star-forming 
regions after a few Myr (Hartmann, Ballesteros-Paredes & 
Bcrgin 2001), and that all clouds of comparable masses 
tend to have comparable levels of turbulence (Heyer & 
Brunt 2004). If the turbulence in clouds were decaying, 
one would expect that clouds of a given mass would ex- 
hibit a large scatter in their turbulence levels, contrary to 
what is observed. The "universal" behavior of the tur- 
bulent level in clouds reported by Heyer & Brunt (2004) 
also suggests that the clouds are part of the global Galac- 
tic turbulent cascade, in which the key process is a sta- 
tistically stationary energy transfer among scales, analo- 
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gous to the classical Kolmogorv cascade. Thus, driven- 
turbulence simulations may be a somewhat better approx- 
imation to real clouds than decaying ones, although the 
standard Fourier-driving scheme is likely not to be the 
best model of the true injection mechanism. Simulations 
with more realistic driving schemes are clearly needed. 

Concerning the subcritical vs. supercritical dichotomy, 
theoretical arguments suggest that as a cloud is becom- 
ing predominantly molecular and self-gravitating, it is 
also becoming supercritical (McKee 1989; Hartmann, 
Ballesteros-Paredes & Bergin 2001). Recent observational 
evidence remains inconclusive, and tends to suggest that 
both kinds are realized, with probably some preponder- 
ance of supercritical ones (Crutcher 1999; Bourke et al. 
2001; Crutcher 2004), although the uncertainties are large. 
The safest assumption at this point appears to be that 
both regimes are realized in molecular clouds. 

Thus, the global picture that emerges is that a distri- 
bution of magnetic field strengths exists in the ensem- 
ble of molecular clouds (Ballesteros-Paredes & Vazquez- 
Semadeni 1997), and that the SFE decreases monotoni- 
cally as the mean field strength in the clouds is larger. This 
appears to be a continuous trend, rather than a sharp di- 
chotomy between sub- and supercritical regimes, as was 
the case in the standard model of star formation (Shu, 
Adams & Lizano 1987; Mouschovias 1991). The role of 
AD would then mostly be to just allow subcritical clouds 
to participate in the star formation process, with the low- 



est SFEs of the spectrum. A confirmation of this picture 
will require a systematic study of the SFE in 3D, driven 
simulations including AD, to be presented elsewhere. 

Our second result is that the number of collapsed objects 
appears to decrease, and their median mass appears to in- 
crease with increasing mean field strength. This result is 
consistent with the observation that in the non-magnetic 
case a single clump appears to form several collapsed ob- 
jects, while in the magnetic cases a clump appears to form 
a single object, as can be seen in the animations presented 
in Paper I. This observation goes in line with the findings 
of Heitsch, Mac Low & Klessen (2001), who noticed that 
non-magnetic simulations with large-scale driving tended 
to form clusters of collapsed objects, while magnetic cases 
tended to form the collapsed objects in a more scattered 
fashion. If collapsing objects in the magnetic case are more 
massive, then clumps of a given mass can form fewer ob- 
jects. Detailed analysis of the evolution of individual cores 
will be necessary to test these possibilities. 

This work has received partial financial support from 
grants CONACYT 36571-E to E.V.-S. J.K. was supported 
by ARCSEC, of the Korea Science and Engineering Foun- 
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(funded by KASI and ARCSEC) and at CRyA (funded by 
the above CONACYT grant). 
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Fig. 1. — Evolution of the accreted mass fraction for the four sets of runs considered. Solid lines denote fi = oo, dotted lines denote /* : 
and dashed lines denote fi = 2.8. 
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Fig. 2. — Masses of the collapsed objects (objects with densities n > 500no) versus the inverse of the simulation mass-to-flux ratio for the 
four sets of runs. The left vertical axis gives the masses as fractions of the total mass in the simulation, and the right axis gives them in solar 
masses, according to the adopted normalization. 
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Fig. 3. — Normalized density histograms for all the simulations immediately before gravity is turned on. The line coding is as in fig. 1. The 
dot-dashed line in the "Paper I" panel denotes the subcritical case with /i = 0.9. 



